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1. SUMMARY 


Numerical and analytical micromechanical models are presented to predict the thermo- 
elastic behavior of a textile composite. In the numerical model, the unit-cell is discretized 
with finite elements, and periodic boundaiy conditions are imposed between opposite faces of 
the unit-cell. For a thin textile composite, stress gradients effects through the thickness are 
demonstrated. The consequent difference in the stiffness and strength behavior of thick and 
thin composites are discussed. The numerical model is implemented to predict the 3-D 
thermo-elastic constants for a thick textile composite, and the plate thermo-mechanical 
properties for a thin textile composite. The numerical model is extended to compute the 
thermal residual microstresses due to processing and to predict the composite failure 
envelopes. An analytical model - Selective Averaging Method (SAM) - is proposed, which is 
based on a judicious combination of stiffness and compliance averaging to estimate the 3-D 
elastic constants. Both the models are tested and verified for several examples by comparing 
the stiffness properties with elasticity solutions and available results. 


2. INTRODUCTION 


The increasing demand for lightweight yet strong and stiff structures has lead to the 
development of fiber reinforced composites. These materials are not only used in the 
aerospace industry but also in variety of commercial applications like automobile, marine and 
biomedical applications. However the manufacture of fibrous laminated composites from 
prepregs is labor intensive. Besides, fibrous laminated composites lack through the thickness 
reinforcement, and hence have poor interlaminar shear strength. Recent developments in 
textile technology shows some promise in overcoming the above limitations. Three- 
dimensional woven and braided composites provide multidirectional reinforcement, thus 
enhancing the strength and stiffness in the thickness direction. Textile manufacturing 
processes such as weaving and braiding in conjunction with resin transfer molding are also 
suitable for the production of intricate structural forms at a reduced turn-round time. 


*Work done on grant at the University of Florida, NAG-1-1226. 
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With the advancements in aforementioned technologies there is a need to develop 
scientific methods of predicting the performance of the composites made using the above 
processes. There are numerous variables involved in textile processes besides the choice of 
the fiber and matrix materials. This, for example, includes the number of filaments in the yarn 
specified by the yam linear density and the yam architecture within the unit-cell determined 
by the type of weaving or braiding processes. Thus there is a need for analytical/numerical 
models to study the effect of these variables on the textile composite behavior. Ideally a 
structural engineer would like to model textile composites as a homogeneous anisotropic 
material - preferably orthotropic - so that the structural computations can be simplified, and 
also the existing computer codes can be used in the design. This would require the prediction 
of the effective (macroscopic) properties of the composites from the constituent material 
(microscopic) characteristics such as fiber and matrix properties, fiber-matrix interface 
characteristics and the fiber architecture. This is possible if we assume that there is a 
representative volume element (RVE) or an unit-cell that repeats its self throughout the 
volume of the composite, which seems to be true in the case of textile composites. Ishikawa 
and Chou (1982; 1983a; 1983b), Yang and Chou (1987), Ma, Yang and Chou (1986), have 
proposed several models for estimating the thermoelastic and mechanical properties of woven 
and braided composites. Yoshino and Ohtsuka (1982), Dasgupta et al. (1990) and Whitcomb 
(1991) analyzed the unit-cell of textile composites using three-dimensional finite elements to 
predict the overall macroscopic behavior of the composites. Their models can be used to 
predict both stiffness and strength properties. In the present paper, we have demonstrated 
numerical models to predict the stiffness and strength behavior of textile composites. An 
approximate analytical method is also described to estimate the stiffness properties of a textile 
composite. 


3. FINITE ELEMENT MODELS FOR THERMO-MECHANICAL PROPERTIES 


In this section, we demonstrate micromechanical models utilizing finite elements to 
predict the effective stiffness properties and coefficients of thermal expansion (CTE's) for a 
textile composite. The macroscopic properties of the composite are determined at a scale 
much larger than the dimensions of the unit-cell, but comparable to the dimensions of the 
structural component. The average stresses at a point at the structural scale will be called the 
macroscopic stresses or macrostresses. The actual stresses at a point at the continuum level 
will be called the microscopic stresses or microstresses .To distinguish the macroscopic 
deformations and stresses from their microscale counterparts - a superscript "M" will be 
used to denote the macroscopic deformations and stresses. 


3.1 Unit-Cell Analysis for Three-Dimensional Elastic Constants 


The unit-cell analysis assumes that the material is subjected to a uniform state of strain 
in the macroscopic sense. The average stresses required to create such a state of strain is 
computed from the finite element model of the unit-cell. In the microscale, all unit-cells have 


626 


I 


identical stress and strain fields. Continuity of stresses across the unit-cell then requires that 
tractions be equal and opposite at corresponding points on opposite faces of the unit-cell. 
Since the displacement gradients are constant for a homogeneous deformation, the 
displacements at corresponding points on opposite faces of the unit-cell differ only by a 
constant. 


Consider a rectangular parallelepiped as the unit-cell of the three-dimensional textile 
composite. The edges of the unit-cell are assumed to be parallel to the coordinate axes x v x 2 
and x 3 , with unit-cells repeating in all three directions. The length of the unit-cell in the x, 
direction is defined as L r On the macroscale the composite is assumed to be homogeneous 
and orthotropic and the composite behavior is characterized by the following constitutive 
relation : 





( 1 ) 


where {a"} and {e w } are the macroscale stresses and strains respectively; {a s } and [C] are 
the macroscale CTE's and orthotropic elasticity matrix to be determined; AT* 1 is a uniform 
temperature difference throughout the unit-cell. 

3.1.1 Periodic Boundary Conditions 

The periodic BC's consist of the periodic displacement boundary conditions which 
ensure the compatibility of displacements on opposite faces of the unit-cell, and the periodic 
traction boundary conditions to enforce the continuity of stresses. A macroscopically 
homogeneous deformation can be represented as 

u* 1 = H..XJ ij =1,2,3 (2) 

where H (j are the displacement gradients. Then the periodic displacement boundary conditions 
to be imposed on the faces x t =0 and x t =L { are 

OyX 2y x 3 ) = H u L t 

= H 12^*2 (3) 
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The traction boundary conditions to be imposed on the faces x~0 and x~L { are 

,) = -F,(x„ O^j) < 4) 

F t (x lt K 2i L 3 ) = -F.(x lt K 2 , 0) 

The above periodic BC's are imposed in the finite element model by using multi-point 
constraint elements or by using transformation equations to eliminate the constrained 
displacements (Cook et al., 1989). Both methods require a finite element model with 
corresponding nodes on opposite faces of the unit-cell. 

3.1.2 Determination of Three-Dimensional Elastic Constants and CTE's 

The unit-cell is discretized with three-dimensional finite elements such that opposite 
faces of the unit-cell have identical nodes. Periodic displacement and traction boundary 
conditions are enforced between the corresponding nodes. The periodic displacement BC's are 
imposed such that only one of the components of the macroscopic strains is non-zero; and the 
uniform temperature difference AT M is set to zero. Then, the difference in displacements 
between corresponding points on opposite faces of the unit-cell will be equal to that in a 
homogenous continuum subject to the same deformation. The average stresses (macroscopic 
stresses) required to create such a deformation are obtained from the finite element results. 
Substituting the macroscopic stresses and strains in the composite constitutive relation Eqn. 

(1) the stiffness coefficients in the column corresponding to the non-zero strain can be 
evaluated. This procedure is repeated for other macroscopic strain components (keeping the 
temperature difference zero) to obtain the entire stiffness matrix [C]. The orthotropic elastic 
constants of the composite material can be easily determined by inverting the stiffness matrix, 
and comparing the compliance coefficients with that of an orthotropic material. 

To compute the six CTE's, a finite temperature change T a is applied to all the elements 
in the unit-cell; and periodic displacement BC's are imposed such that all the macroscopic 
strain components are zero. Then the composite constitutive relation Eqn. (1) will reduce to 

{o M } = -[C]{<x*} T 0 (5) 

The macroscopic stresses for such a deformation are computed as described below. Knowing 
the stiffness matrix [C], the composite CTE’s are found as 

M - -jric r l {°"} <«> 

1 o 

Table 1 presents the non-zero displacement BC’s imposed on the unit-cell to obtain [C] and 
the CTE's {a 1 }. 

The macroscale stresses for a given deformation state can be found by one of the 
following two methods. In the first method, the macroscale stresses are obtained by averaging 
the nodal forces on each face of the unit cell. For example, the macroscale o,, M can be 
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obtained as 


o"= (7) 

r 2^3 

where F/ nJ is the nodal force in the x, direction at the nth node, and Z„ denotes summation 
over all nodes on the face x,=L,. Alternately, the macroscopic stresses can be computed by 
volume-averaging the corresponding microstress component in the unit-cell. Then the 
macroscale (S U M is obtained as 

ofi = -pj* o u Qc&*)dV (8) 

v v 

where V is the volume of the unit-cell. The microstresses are computed at the quadrature 
points, and numerically integrated over the volume in each element of the unit-cell. 

3.1.3 Results and Discussion for 3-D Elastic Constants and CTEs 

The above procedure was demonstrated for the following materials: 

Example 1. isotropic material 

Example 2. bimaterial medium - both materials are assumed isotropic 
Example 3. unidirectional composite with identical poisson ratios for fiber and matrix - 
fiber and matrix materials are isotropic 

Example 4. unidirectional composite with different poisson ratios for fiber and matrix - 
fiber and matrix materials are isotropic 

Example 5. plain-weave textile composite (yam geometry and properties obtained from 
Dasgupta et al., 1990) 

Example 6. plain-weave textile composite (yam geometry and properties obtained from 
Naik, 1994) 

Example 7. 5-hamess satin weave (yam geometry and properties obtained from Naik, 1994) 
For the textile composite examples, i.e, examples 5-7, the yam is assumed to be transversely 
isotropic and the matrix material is assumed isotropic. The constituent material properties for 
the examples are listed in Table 2. 

A 3-D finite element code called /uTE%-20 (pronounced as mutech ) was written and 
implemented for the seven examples to compute the homogeneous elastic constants and 
CTE's. The unit-cell was divided into uniform eight noded hexahedral elements. The element 
stiffness matrix (K e ) was formulated as 

K e = f B t CB dV e 

Jy > 

+ 1 + 1+1 

= fff B t CB\J\ dld^di (9) 

- 1 - 1-1 

■LEE 

i=l j = 1 *=1 

where V" is the domain of the element, B is the strain-displacement transformation matrix, C 


629 


is the elasticity matrix, N is the number of Gauss points used for integration, W is the 
Gaussian integration weight factor and |/| is the determinant of the Jacobian. The material 
property at each Gauss point (£,T|,0 was determined, and the corresponding elasticity matrix 
was used to perform the volume integration. The element stiffness matrix thus represented the 
averaged properties of the constituent materials in that element. The computed elastic 
constants for the seven examples are listed in Tables 3-5. 


The bimaterial medium consisted of two different layers of equal thickness in the xy- 
plane alternating in the z-direction (Fig. 1). The effective Youngs moduli, Poisson ratios and 
CTE’s of the bimaterial medium were derived exactly as described below. The constitutive 
relation (considering only the normal stresses) for each layer was defined as 
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( 10 ) 


where the superscript refers to the layer number. To derive C n , C 12 and C ]3 for the bimaterial 
medium, £ XX M was assumed as one; and e yy M and e u M as zero.The assumption of =0 
implied that £ u ‘= -ej. The following constraints were, in addition, imposed across the 
bimaterial interface : 


1 2 _ M 

€ x* e xc e x* 


1 = c 2 = e M 

e yy e yy e yy 


'zz 


2 u 

a zz ~ a zz 


( 11 ) 


From the above interfacial constraints and Eqn. (10), the stresses in each layer were 
computed. The stresses in the layers were volume-averaged to yield the corresponding 
macroscopic stresses, i.e, G„ M , c yy M and O zz M . Since £ XX M was equal to unity, the computed 
macroscopic stresses were identical to the stiffness coefficients C I2 , C I2 and C I3 . A similar 
procedure was followed to find the remaining stiffness coefficients and CTE's for the 
bimaterial medium. The inplane shear modulus of the bimaterial medium was computed as 
G x = (G,+ G 2 )/2 knowing that the shear strain was uniform in both layers. The isostress 
assumption was used to derive the transverse shear modulus as G xl - 2G,G 2 / (G,+G 2 ). It was 
found that /uTEx-2 0 results were identical to the elasticity results for the bimaterial medium 
(Table 3). 

Table 4 presents the elastic constants and CTE’s for the two unidirectional composite 
examples. The unidirectional composite unit-cell is shown in Fig. 2. The unidirectional 
composite properties were compared with available analytical solutions. The rule of mixtures 
formulae were used to predict E L and v Lr ; the Halpin-Tsai equations (Halpin and Tsai, 1969) 
for E v G lt and and Schapery's expressions (Agarwal and Broutman, 1990) for the thermal 
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coefficients a L and a T . The relations for E L , \ LT and a L are exact when the poisson ratios are 
identical for the fiber and the matrix. The reason for the discrepancy in example 3 is that a 
point-location subroutine in /uTE%-20 identified a fiber volume fraction of 0.595, whereas the 
actual fiber volume fraction was 0.6. 

The elastic constants for the three textile composite examples are presented in Table 5. 
Figures 3 & 4 illustrate the weave patterns for the plain-weave (examples 5 & 6) and the 5- 
hamess satin weave (example 7) respectively. The properties for example 5 were compared 
with Dasgupta's results for an overall fiber volume fraction (V f ) of 0.26. The yam properties 
were not specified in Dasgupta et al. (1990). So the rule of mixtures and Schapery's 
expressions were used to obtain the yam properties from the given fiber and matrix 
properties. The pTE %- 20 results for examples 6 and 7 were compared with TEXCAD - an 
approximate analytical method developed by Naik (1994). In both the examples the overall 
fiber volume fraction was 0.64. It must be noted that nTE %- 20 will marginally under-predict 
the stiffness moduli - since the yam cross-section in the numerical model is approximated as 
a polygon inscribed within the actual cross-sectional area. Consequently the yam/fiber volume 
fraction in the numerical model will always be lesser that the theoretical volume fraction. 


3.2 Stress Gradient Effects 


The methods explained in Section 3.1 assume that the unit-cells exist in all the three 
directions. This will be tme in the case of thick textile composites. However there are many 
applications in which thin composites are used. In fact in order to take advantage of the 
properties of composites, the structures have to be made of thin plate like members with 
stiffeners for load transfer. In such cases there will be fewer unit-cells in the thickness 
direction. Thus the free surface effects will be predominant. There will be severe stress 
gradients through the thickness, and they will have an influence on the apparent stiffness and 
strength of the structure. 

The following simple example will illustrate the stress gradient effects on stiffness. 
Consider a layered medium consisting of alternating layers of materials of equal thickness 
with Young's moduli E, and E 2 respectively (Fig. 5a). Any micromechanical model would 
predict that the medium can be considered as a homogeneous orthotropic material at 
macroscale and also the effective Young's modulus in the longitudinal direction is (E,+E 2 )/2, 
and there would not be any bending-stretching coupling in the principal material direction. 
However, if we consider a bimaterial beam consisting of the same two materials (Fig. 5b), we 
will find that there is a bending-stretching coupling, and also the flexural rigidity cannot be 
predicted from the Young's modulus of the homogeneous orthotropic medium and the total 
beam thickness. The bimaterial beam has properties and behavior different from the 
corresponding infinite medium. This phenomenon is observed in the transverse shear behavior 
also (Sankar and Marrey, 1993). A similar behavior is also expected in thin textile composites 
where there are fewer unit-cells in the thickness directions, and the unit-cells are not 
subjected to a macroscopically homogeneous state of deformation as assumed in Section 3.1. 
The effect of stress gradients on stiffness and strength of thin textile composites are discussed 
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in Marrey and Sankar (1993a). 

One method of overcoming the above difficulties in thin textile composites is to model 
the composite as a plate/beam, and compute the structural stiffness properties (eg., [A], [B] 
and [D] of the plate) directly from the unit-cell analysis instead of the continuum stiffness 
properties (Young's modulus, shear modulus etc.). In the following sub-sections we illustrate 
these concepts - first for a thin textile composite modeled as a beam, and then for a textile 
composite plate. The purpose of the beam model is to present the issues involved in 
computing the structural stiffness coefficients. Further the periodic BC’s are different from 
those in the continuum model. 


3.3 Unit-Cell Analysis for Beam Thermo-Mechanical Coefficients 


The textile composite beam is assumed to be in the xz-plane with unit-cells repeating 
in the jc-direction. A state of plane strain parallel to the xz-plane is assumed. On the 
macroscale it is assumed that the beam is homogeneous and its behavior can be characterized 
by the following beam constitutive relation : 
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where P, M and V are the axial force, bending moment and transverse shear force resultants 
respectively; [AT] is the symmetric matrix of beam stiffness coefficients, Eg , k** and Yo are 
the midplane axial strain, curvature and transverse shear strain respectively; a P , a v and a M 
respectively are the thermal expansion, thermal shear and thermal bending coefficients. The 
midplane deformations are related to the midplane axial displacement u 0 , transverse 
displacement w, and rotation \|/ as: 


*"a? Yo * a. 


( 13 ) 


Actually K u , K n , K 22 and K n are similar to the laminate stiffness coefficients A u , B u , D n and 
y?A 55 respectively. The beam constitutive relation in Eqn. (12) can also be expressed in terms 
of compliance coefficients as 
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As discussed earlier, the unit-cell analysis assumes that all the unit-cells are subjected 
to identical stress and strain fields, for a given state of loading. This is true in the case of 
constant axial force (P) and constant bending moment (M) in the beam. However, when a 
shear force (VO is applied to the beam, the shear force will give rise to the building up of 
bending moment at every cross-section, such that V = - (dM/dx). This situation created 
difficulties in estimating the shear stiffness of the beam accurately. 

The detailed procedures for evaluating [A"] and the beam CTE's are explained in 
Sankar and Marrey (1993) and Marrey and Sankar (1993b). However the principles involved 
in finding the beam stiffness matrix [AT] are described briefly for the sake of completion. 

Three linearly independent deformations are applied to the unit-cell namely. 

Case (i) unit axial strain (e 0 = 1, K = 0, y 0 = 0) 

Case (ii) unit curvature accompanied by transverse deflection such that the transverse 

shear strain was zero (e 0 = 0, k = 1, y 0 = 0) 

Case (iii) unit transverse shear strain (e 0 = 0, k = 0, y 0 = 1 ) 

The periodic displacement boundary conditions for the three unit deformations are given in 
Table 6. The temperature change, AT, is assumed to be zero. The top and bottom surfaces of 
the unit-cell are considered as free surfaces. For each case, the axial force P, the bending 
moment at the center of the unit cell M, and the shear force V are computed. Since the 
bending moment varies linearly along the unit cell, the bending moment at the center will be 
the average of the bending moments at the left and right ends of the unit cell. By substituting 
the values of P, M and V in Eqn. (12), one can evaluate the stiffness coefficients. 

3.3.1 Beam Results and Discussion 

The procedures to obtain the beam coefficients were implemented for the following 

cases: 

(a) an isotropic beam; (b) a bimaterial beam with isotropic layers of equal thickness; (c) a 
plain weave textile composite beam where the yarn is assumed to be transversely isotropic 
and the matrix is isotropic. The properties of the constituent materials for all the cases are 
listed in Table 7. The dimensions of the unit-cell and the yam architecture for the textile 
beam were taken from Yoshino and Ohtsuka (1982). The same unit-cell dimensions (length of 
3.6 mm and height 1.8 mm) were also used for the isotropic and bimaterial cases. 

The unit-cell of the beam was discretized using eight-noded isoparametric plane strain 
finite elements. The finite element mesh for the isotropic unit-cell and the plain weave unit- 
cell were identical except that different material properties were used. The deformed plain 
weave unit-cell under various independent loading conditions is shown in Fig. 6. 

The stiffness and thermal coefficients for the three beams examples are shown in 
Table 8. The results for the isotropic and bimaterial beams were compared to exact beam 
theory solutions. Exact shear correction factors - 0.833 for the isotropic beam and 0.555 for 
bimaterial beam were used in the beam theory solution to compute the shear stiffness 
(Whitney, 1973). It can be seen that the beam unit-cell analysis is able to predict the axial 
and bending stiffness coefficients {K n and K 22 ) very accurately. As expected the shear 
stiffness (K n or A 55 ) predictions have errors, but they are very minimal. The textile beam 
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stiffness coefficients were also estimated using a procedure similar to the mosaic model 
(Ishikawa and Chou, 1982). They are compared with the coefficients obtained from the unit- 
cell analysis as shown in Table 8. It can be seen that the mosaic model predicts K 33 
reasonably well. The reason for the lack of agreement in K n and K 22 can be attributed to the 
fact that a major portion of the yam is modeled as a 0° laminate in the mosaic model, which 
tends to over-predict the axial and flexural stiffness. 

The shear modulus of the plain-weave beam G xz was found by imposing periodic 
boundary conditions between the top and bottom surfaces, and left and right ends of the unit- 
cell. This would yield the apparent shear stiffness as G xz h= 5.53 x 10 6 Nm' 1 - whereas the 
actual shear stiffness is 9.21 x 10 6 Nm' 1 (K 33 in Table 8). The Young’s modulus of the textile 
beam E x may be extracted from K,„ as K,/h, which would yield E= 15.42 GPa. If this value 
of E z were used to predict the homogeneous flexural stiffness as D n =E x h 3 /12, we would 
obtain D,, as 7.50 Nm - whereas the actual flexural stiffness is equal to 5.41 Nm. The same 
idea holds for the beam thermal coefficients also. The beam CTE's ,i.e, a P , a M and a v cannot 
be predicted from the corresponding continuum CTE's. Table 9 shows the disagreement for 
the plain-weave example between the beam CTE s obtained directly , and the beam CTE s 
predicted from the corresponding continuum CTE's. It may be noted that the continuum 
model would always predict the thermal expansion coefficient cl p as ( X x , and the thermal 
bending coefficient a M as zero. This underscores the importance of the present analysis for 
predicting the beam stiffness properties for a thin textile composite directly. 


3.4 Unit-Cell Analysis for Plate Thermo-Mechanical Coefficients 


In this section we describe a procedure to find the stiffness and thermal properties of a 
textile fabric modeled as a structural composite plate. The textile composite plate is assumed 
to be in the xy-plane with unit-cells repeating in the a and v directions. The lengths of the 
unit-cell in the x- and y-directions are assumed to be a and b respectively and the unit-cell 
thickness as h. On the macroscale the plate is assumed to be homogeneous and the plate 
behavior is characterized by the plate constitutive relation: 
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where e i0 M , y i0 M and k, m are the midplane axial strain, shear strain and curvature; a/’ and 
are the plate thermal expansion and bending coefficients; N t and M, are the axial force and 
bending moment resultants respectively in the homogeneous plate. The midplane strains and 
curvatures are related to the midplane displacements and rotations as: 


M _ dUp M _ U _ fop 

ax' ay' ay 

M _ M _ 

~ dy' ay 


dx 

+ d% 
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(16) 


(17) 


3.4.1 Unit-Cell Boundary Conditions 


The plate uicrmo-ir.cchar.ical properties were obtained hv modeling the unit-cell with 
eight-noded brick elements and subjecting the unit-cell to linearly independent deformations. 
The unit-cell was subjected to minimum support constraints to prevent rigid body rotation and 
translation. The top and bottom surfaces of the plate were assumed to be free of tractions. 

The faces x=0 and x=a had identical nodes in the finite element model, and so did the pair of 
faces y-0 and y=b. The identical nodes on opposite faces of the unit cell were constrained to 
enforce the periodic BC's. The traction boundary conditions on the lateral faces of the unit- 
cell were: 


FfWd = -F.( Oyj), FfxM = -FfrA z), i= W ( 18 > 


The periodic displacement BC's enforced for the deformations are presented in Table 10. 
3.4.2 Determination of Plate Stiffness Coefficients and CTE's 


Linearly independent deformations are applied to the unit-cell such that only one of 
the six components of deformation is non-zero (first six cases in Table 10). The temperature 
difference is set to zero for all six cases. It must be noted that the applied deformations must 
ensure that the transverse shear strains, y xz M and y n M are zero where 
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The force and moment resultants can be obtained by averaging the nodal forces on each face 
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of the unit cell. For example, on the face x=a the force and moment resultants are computed 
using the relations: 


N x = EF®(flj,z), 
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where F/ ,J and F y w are the nodal forces in the x and y directions at the i' h node and 'n' is the 
total number of nodes on the face. The force and moment resultants can also be computed by 
averaging the microstresses over the unit-cell volume. Then the resultants on the face x=a are 
obtained as 


N x = N xy = 

M r, = j- b fzoJxj*)dV = -Lfz^ixjJdV 


( 22 ) 
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Substituting the values of the deformation and the force resultants in the plate constitutive 
relation, Eqn. (15), the stiffness coefficients in the column corresponding to the non-zero 
deformation can be computed. This procedure is repeated for other deformation components 
to obtain all the stiffness coefficients. 


To predict the CTE’s, the plate unit cell is subject to a uniform temperature difference, 
given by AT = T 0 . In the finite-element model, periodic displacement BC's are applied such 
that all six components of the deformation are zero (seventh case in Table 10). The force and 
moment resultants are computed using one of the procedures described above. The thermal 
expansion coefficients a p , and thermal bending coefficients are then obtained from the 
relation: 
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3.4.3 Results for Plate Stiffness Coefficients 


The plate [A], [fi], [D] matrices and CTE’s were found for the seven examples listed in 
Table 2 by implementing the finite element code /uTEx-2 0. The plate properties for the 
isotropic and bimaterial cases are presented in Tables 11 and 12 respectively. The bimaterial 
plate properties were also computed using the lamination theory for two plies, and from the 
3-D elastic constants (Table 12). For example the coefficient D n is obtained from the 3-D 
elastic constants as 
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n 
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The finite element results for the bimaterial case were exact, i.e, identical to the results 
obtained with the two-ply lamination theory. The [A] and [£>] matrices computed from the 
bimaterial 3-D constants were found to be in good agreement with the two-ply lamination 
theory only because both the layers were equal in thickness. In-general, however the plate 
properties obtained from the 3-D elastic constants would be different from the two-ply 
lamination theory results. The plate properties for the unidirectional composite examples are 
presented in Table 13 and for the textile examples in Table 14. In all the examples it was 
found that the plate properties, especially [5], D u , {a p } and {p''} could not be predicted 
from the corresponding 3-D elastic constants. 


4. FINITE ELEMENT MODELS FOR STRENGTH PROPERTIES 


In the previous sections, we had demonstrated methods to model a general textile 
composite either as a three-dimensional material nr as a thin nlate/beam to Dredict their 
corresponding thermo-mechanical coefficients. In this section we extend the same numerical 
models to compute the thermal residual stresses due to processing in the yams and the matrix. 
Then the numerical models are used to study the strength behavior of the composite 
by predicting the failure envelopes for thin and thick textile composites. 


4.1 Thermally Induced Residual Microstresses 


The mismatch in the CTE's for the constituent materials in the composite induces the 
residual microstresses in the yam and matrix. The difference between composite curing 
temperature and room temperature then serves as the driving force to create these 
microstresses.The microstresses in the vicinity of the yam-matrix interface are particularly 
important as they could lead to failure due to debonding. Since composites designed for high 
temperature applications are fabricated at higher temperatures, the residual microstresses 
become particularly relevant in the strength considerations of such composites. 

4.1.1 Determination of Residual Microstresses 

Let T a be the difference between room temperature and the composite fabrication 
temperature. Since the composite is stress free at the fabrication temperature, which is above 
the room temperature, T a is generally negative. The residual microstresses in the yam and the 
matrix are obtained by superposing the microstresses due to the two load cases as explained 
below. In the first load case, the unit-cell is constrained from expanding by fixing the comer 
nodes of the unit-cell and enforcing zero displacement difference between corresponding 
nodes on opposite faces of the unit-cell (periodic displacement BC's). A temperature 
difference T a is applied to all elements in the finite element model. This is exactly the same 
problem we solve for finding the three-dimensional CTE's. The applied boundary conditions 
mean that all the macroscopic strain components are equal to zero ( {t^ 1 } = 0 , A I* 1 = T 0 ). 
Then the corresponding macroscopic stresses required to restrain the unit-cell expansion are 
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given by 


{«") - ~[C]{a*}T 0 1261 

In the second load case, deformations are applied so as to reverse the macroscopic stresses 
developed in the first load case. This can be accomplished by imposing the deformations 
{e M } = {a'} T 0 and AT M =0. It can be noted that the macroscopic stresses developed in the 
second loading case [C]{a'}r o are equal and opposite to the macrostresses in Eqn. (26). The 
microstresses from both load cases are superposed to obtain the residual stresses due to free 
thermal expansion. 

The same idea can be extended to finding the residual microstresses in the plate 
model. Then the deformations to be applied in the first load case are {e M } =0, { } = 0 and 
A7" = T 0 ; and the deformations in the second load case are {e"} = {a p } T a , {k"} = {J^} 
T o and AT" = 0. The residual microstresses were computed for the plain-weave textile beam 
at the Gaussian center of the elements in the unit-cell. Figure 7 shows the thermal stress 
contours for a u and X xz . The composite curing temperature was assumed to be 150°C 
above room temperature. 


4.2 Strength Modeling of Textile Composites 


There are many failure criteria or strength theories for unidirectional fiber composites. 
This for example includes (Agarwal and Broutman, 1990) maximum stress theory, maximum 
strain theory, Tsai-Hill theory. Even though failure of a material is a very complex 
phenomenon, engineering strength theories such as mentioned above have been found to be 
useful in design. The interpretation of strength values obtained from such theories are 
different for different materials. For example in metal matrix composites the failure envelope 
obtained using the above theories will correspond to the initial yield surface (Dvorak et al., 
1973). In graphite/epoxy composites the failure theories can be used to predict fiber or matrix 
failure. In the present study our intent is to explore the possibility of developing such failure 
criteria for textile composites. 

4.2.1 Determination of Composite Failure Envelope 

Our approach is similar to that used by Dvorak et al. (1973). A state of homogeneous 
deformation, corresponding to each of the six macrostrain components, are independently 
applied to the unit-cell by imposing the boundary conditions explained in Section 3.1. For 
each case, the various stress components are computed in the elements in the unit-cell - 
typically at the element Gaussian integration points. These stresses will be referred to as 
microstresses. Assuming linear elastic behavior, the microstresses can be computed for any 
arbitrary combination of the macrostrain components. Since we know the macroscopic elastic 
constants, we can find a relation between microstresses at a point and any arbitrary state of 
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macrostress as: 


{a} = [F] {o w } 


(27) 


[F] can be considered as a matrix of influence coefficients, which is evaluated at the 
integration points of all the elements in the unit-cell. We also assume that the failure behavior 
of the matrix material and the yam is known. For instance, let the failure criterion of the 
matrix be given by [//] { or } matrix = 1. Then the failure criterion for the composite is obtained 
from Eqn. (27) as [//] [F] {o"} = 1. The same idea also applies for the yam. The textile 
composite is assumed to have failed, if there is failure on the microscale in any one of the 
constituent materials - either matrix or the yam. By varying the macrostresses using a 
numerical simulation, failure envelopes can be obtained for the idealized homogeneous 
material. It might be noted that Eqn. (27) can be modified to include the thermal residual 
stress field in the unit-cell as 


(o| = \F\ (28) 

\ -j 

where {a r } is the matrix of thermal microstresses computed at the element integration point 
for a temperature difference of T a . 

4.2.2 Effect of Stress Gradients on Strength 

The strength analysis for a three-dimensional composite can extended for thin 
composites using the plate model. As mentioned in Section 3.2, the macroscopic stresses will 
not be homogeneous thorough the thickness in such composites. Then the composite failure 
will be determined by the stress gradients through the thickness, represented by the averaged 
force resultants (F) and moment resultants ( M ). The composite failure criterion will be of the 
form 


[H][F\ 



= 1 


(29) 


Thus the failure envelope of the composite will be in the six-dimensional space of the force 
resultants and the moment resultants. The above procedure was demonstrated using the beam 
model, for the case of a plain-weave textile composite. 

4.2.3 Beam Failure Envelope Results 

The strength properties used for the constituent materials in the beam are as follows: 
Yam: a/=1725 MPa, o L c =1366 MPa, o/=42 MPa, a r c =230 MPa, x L1 =95 MPa 

Matrix: o r = 70 MPa, a c =100 MPa 

where the superscripts T and 'C' refer to the tensile and compressive strengths respectively. 
The maximum principal stress criterion was used to determine the matrix failure and the 
maximum strain theory for unidirectional fiber composite was used for the yam. Both 
continuum and structural failure envelopes were developed. Figures 8 and 9 depict the 
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structural failure envelopes in the P-M space. If we assume that the beam is made of a 
homogeneous but orthotropic material with properties as predicted by the continuum model, 
then one can derive structural failure envelopes from the continuum failure envelopes using 
simple beam theories. The derived structural failure envelopes are compared with that 
obtained from direct micromechanical analyses in Figs. 10 and 11 respectively. One can note 
that the continuum failure criteria are very conservative for the case of a thin beam. 


5. ANALYTICAL MODELS FOR STIFFNESS PROPERTIES 


The complex yam architectures in a textile composite make numerical modeling of 
the unit-cell extremely difficult. Analytical methods are popular because they are easy to 
model and suitable for performing parametric analysis. However these methods are 
approximate because they assume a certain state of stress or strain in the unit-cell. Averaging 
the stiffness or compliances of the matrix and the inclusion has been used to estimate the 
bounds of effective elastic properties of the composite. Essentially the stiffness averaging 
assumes a state of uniform strain in the composite (isostrain), and compliance averaging 
assumes a state of uniform stress (isostress) in the matrix and inclusion. In fact the rule of 
mixtures expressions for estimating the effective properties of a unidirectional composite is 
based on such averaging schemes. These methods are too simplistic, because the state of 
stress/strain in a textile composite under a uniform macrostress is much more complex. We 
propose a scheme of selective averaging in which both stiffness and compliance coefficients 
can be averaged selectively depending on a more realistic assumption of either isostress or 
isostrain. 


5.1 Selective Averaging Method (SAM) 


Consider a rectangular parallepiped of dimensions axbxc as the unit-cell. To find the 
first column of the effective stiffness matrix, the unit-cell is divided into slices (mesolevel) of 
thickness d parallel to the yz-plane (Fig. 12a). Each slice is further sub-divided into elements 
(microlevel) as shown in Figs. 12(b) and 12(c). In this section, to distinguish between the 
macrolevel, mesolevel and microlevel properties, an over-tilde is used to denote the mesolevel 
properties, and a superscript "M" is used to denote the macrolevel properties. For example, 
[C M ], [C] and [C] will represent the macrolevel, mesolevel and microlevel stiffness 
respectively. The unit-cell is subjected to a deformation such that all macrostrains except z xx M 
are equal to zero and 1 . It is assumed that the mesolevel and microlevel strains, 

corresponding to the zero macrostrains, are negligible. 


The average stiffness of a slice can be obtained based on the isostrain assumption within the 
slice ( ejx.y.z) = (*) ) as : 
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(31) 


C n (x)=j-ff C n {x,yj)dydz 

where C u (x,y,z) is the element stiffness coefficient transformed to the unit-cell coordinates. 
The stiffness of the slices are averaged on the macrolevel based on the isostress assumption, 
i.e., <5^ (x) = o xx M . Thus the first column of the effective stiffness matrix can be computed 
using the following relations: 


1 _1 r 1 

c" «■'£„« 


dx 



1 

abc 



C l7 (w) dxdydz 


(i=2,...,6) 


(32)) 


A similar procedure can be implemented to determine the second and third columns of the 
homogenous stiffness matrix [C M ]. 

A slightly different averaging scheme is used when the unit-cell is subjected to shear 
strains on the macrolevel. Consider the case where the unit-cell is subjected to unit y n M at 
macroscale. Assume that all the other components of strain at the macrolevel, mesolevel and 
microlevel are zero, i.e., 


M 


= e. = e. 


= 0 


i* 4 


(33) 


where e 4 = J yz . We also assume that the shear stress is constant in a slice such that X yz (x,y,z) 
= t yz (x). The shear compliance of a slice can then be obtained by averaging the shear 
compliances of all the elements in the slice as: 


1 

CJx) 



dy dz 

C u (xyj.) 


(34) 


The fourth column of the stiffness matrix Cj 4 is obtained under the assumption that all the 
slices are under a state of constant shear strain: 




c u (w) 


C i4 (x,yj) dxdydz 


(i 1,...,6) 


(35) 


A similar procedure is used to determine the fifth and sixth columns of [C w ]. 
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5.2 SAM Results 


A code was written in FORTRAN to implement SAM. The code (juTE%-\0) was 
used to predict the elastic constants for the seven examples, whose constituent material 
properties are listed in Table 2. Input to the code were the unit-cell dimensions, yam/fiber 
geometry, constituent material properties, and the number of divisions required to discretize 
the unit-cell in the x, y and z directions. The element stiffness matrix [C] was computed by 
finding the constitutive stiffness matrix for the point at the center of the element, and 
transforming it to the unit-cell coordinate system. 

The results for the bimaterial medium (example 2) are given in Table 15. The 
bimaterial medium consisted of two layers of isotropic materials stacked alternately in the z- 
direction (Fig. 1). It can be observed that SAM marginally under-predicts the bimaterial 
longitudinal and transverse Youngs moduli, while the inplane and transverse shear moduli are 
exact. Table 16 presents the SAM results for two cases of unidirectional composite (examples 
3 and 4). The fiber and matrix had identical poisson ratios in example 3, and different 
poisson ratios in example 4. The elastic constants from SAM were compared with the finite 
element results from Section 3.1 and with analytical solutions for a unidirectional composite 
(the rule of mixtures formulae and the Halpin-Tsai equations). Table 17 compares the SAM 
properties for three textile composites (examples 5, 6 and 7) with the previously computed 
finite element results and other available results. In all three cases the elastic constants 
obtained by implementing SAM were in good agreement with the available results. 


6. CONCLUSIONS 


Micromechanical models were demonstrated to predict the stiffness and strength 
behavior of textile composites. In thin textile composites, the stress gradients through the 
thickness were significant, and hence the composite was modeled as a homogeneous 
plate/beam rather than homogeneous continua. The beam model was first discussed to present 
the issues involved in computing the structural stiffness and strength properties. Then the 
plate model was presented for computing the plate stiffness and plate CTE's, and to predict 
the failure envelope for a thin textile composite. The failure envelope for a thin textile 
composite was described in the space of the force and moment resultants instead of the space 
of the macroscopic stresses. It was shown for various examples, that the plate properties could 
not be predicted from the corresponding continuum properties. 

In the models using finite elements (Section 3), the periodic BC's were enforced for 
continuity of tractions and displacements across the unit-cell boundaries. Therefore the finite 
element results for the stiffness properties could be expected to be very accurate. However 
due to difficulties in mesh generation, the unit-cell for the 3-D continuum and the plate 
models were discretized into inhomogeneous finite elements (elements with one or more 
constituent materials). Thus the stresses in the vicinity of the yam-matrix interface for the 
above models may be approximate. For the beam model, however, the unit-cell was 
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discretized into homogeneous elements. 


The Selective Averaging Method (SAM) discussed in Section 6 was based on a 
combination of isostress and isostrain assumptions. The method was fast and easy to 
implement and suitable for parametric studies such as yam preform selection. Extension of 
SAM to predict the continuum CTE's, plate stiffness and plate CTE's is straightforward and is 
underway. 
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Table 1. Non-zero BC's Imposed to Obtain 3-D Thermo-Elastic Constants 


stiffness coefficients to be obtained 

non-zero EC’s 

first column of [C] (e ; , M = 1 ) 

u } (Lj, x 2 , x 3 ) - Uj(0 y x 2 , x$) — Lj 

second column of [C] ( e 22 m = I ) 

U^Xjy L 2 , X 3 ) - U 2 (Xj, 0y x$) — L 2 

third column of [C] (e 33 M = 1 ) 

U^(x Jy X 2 y Lj) - Uj(Xjy X 2 y 0) — Lj 

fourth column of [C] (y 23 M = 1) 

u 2 (Xjy x 2 , Lj) - u 2 (Xjy x 2 y 0) — Lj/2 
UjfXj, L 2 y Xj) - Uj(Xjy 0 y Xj) — L-/2 

fifth column of [C] (y I3 M = T) 

Uj(Xjy x 2 y Lj) - Uj(xjy x^ 0) — Lj/2 
Uj(Lj t x 2 y Xj) - Uj(0 f x 2 y Xj) — L/2 

sixth column of [C] ( y I2 M = 1) 

ufx,, L 2 , x 3 ) - ufxj, 0, x 3 ) = L/ 2 
u 2 (L } , x 2 , x 3 ) - u 2 (0, x 2 , x 3 ) — L/2 

CTE’s (AT* 1 = 7) 

> 

•s 

II 
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Table 2. Properties of Constituent Materials for Examples 1-7 


Example 1 


E = 10 GPa, v=0.3, a =10x10* /°C 


unit-cell size: 

0.500x0.500x0.256 mm 

Example 2 

layer 1 (E-glass): 

E,=70 GPa, V ,=0.200, a,= 5x10* /°C, V,=0.5 


layer 2 (epoxy): 

E 2 =3.50 GPa, V 2 =0.3 50, a 2 = 60x10* /°C, V 2 =0.5 


unit-cell size: 

0.500x0.500x0.256 mm 

Example 3 

fiber: 

Ej=100 GPa, Vj=0.300, 09 = 10x10* /°C, V j=0.6 


matrix: 

E m =10 GPa, V m =0.300, a n = 100x10* /°C 


unit-cell size: 

10x10x10 fjm 

Example 4 

fiber (E-glass): 

E r 70 GPa, Vj=0.200, Ctj= 5x10* /°C, V j=0.6 


matrix (epoxy): 

E m =3.50 GPa, V m =0.350, a m = 60x10* / °C 


unit-cell size: 

10x10x10 fum 

Example 5 

yarn properties (glass-epoxy): 


E l = 58.61 GPa, E T = 14.49 GPa , G LT =5.38 GPa, \ LT =0.250 


v^O.247, a L =6. 15x10 

*/°C, 0^=22.64x10* /°C, Vj=0.26 1 


matrix properties (epoxy): 


E=3.45 GPa, v=0.37, 

a =69x10* /°C 


unit-cell size: 

1.680x 1 .680x0.228 mm 

Examples 6, 7 

yarn properties (graphite-epoxy): 


E L = 144.80 GPa, E T =11.73 GPa, G LT =5.52 GPa, V LT =0.230 


Vrr=0.300, a L = -0.324x1 O' 6 /°C, 0^=14.00x10* /°C, Vj=0.64 


matrix properties (epoxy): 


E=3.45 GPa, v=0.35, 

a =40x10* /°C 


unit-cell size: 

2.822x2.822x0.2557 mm (Example 6) 



7.055x7.055x0.2557 mm (Example 7) 
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Table 3. 3-D Properties for Examples 1 and 2 using Finite Elements 



Table 4. 3-D Properties for Examples 3 and 4 using Finite Elements 
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Table 5. 3-D Properties for Examples 5, 6 and 7 using Finite Elements 




E X ,E, 

(GPa) 

E z 

(GPa) 

G «> 

G * 

(GPa) 

G xy 

(GPa) 



a/,a/x 

io‘/°c 

a/x 

10‘/®C 

Example 5 
(plain- 
weave) 

fiTE%-20 

(FEA) 

11.55 

6.26 

1.94 

2.12 

0.399 

0.186 

29.75 

75.82 

Dasgupta 

results 

14.38 

6.25 

1.94 

3.94 

0.463 

0.167 

22.5 

86 

Example 6 
(plain- 
weave) 

ftTEx-20 

(FEA) 

58.27 

10.92 

4.32 

4.84 

0.363 

0.097 

1.47 

23.15 

TEXCAD 

64.38 

11.49 

5.64 

4.87 

0.396 

0.027 

1.33 

20.71 

Example 7 
(5-harness 
weave) 

ftTEx- 20 
(FEA) 

62.51 

11.02 

4.45 

4.80 

0.349 

0.047 

1.58 

23.18 

TEXCAD 

66.33 

11.51 

4.93 

4.89 

0.342 

0.034 

1.46 

21.24 


647 
























































Table 6. Periodic Displacement BC's for Beam Stiffness Coefficients 



Table 7. Constituent Material Properties for Beam Examples 


isotropic 

beam 

E = 10 GPa, v = 0.30, a = lOxlO' 6 /°C 

bimaterial 

beam 

Ej = 70 GPa, v ; = 0.33, a, = 23xl0' 6 /°C 
E 2 = 3.5 GPa, v 2 = 0.35, a 2 = 60x1 0‘ 6 /°C 

plain-weave 
textile beam 

yarn: 

Ej = 159 GPa, E 2 = 10.9 GPa, G n = 6.4 GPa, v I2 = 0.38, 

v 2J = 0.38, a, = 0.045xl0' 6 /°C, a 2 = 20.2xl0' 6 /°C where the yam 

direction is parallel to the 7 -axis and 25-plane is the plane of 

isotropy. 

matrix: 

E m = 3.5 GPa, v m = 0.35, a m = 60xl0‘ 6 7°C. 
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Table 8. Comparison of Beam Stiffness Coefficients and CTE's (SI Units) 





k 12 

K 22 

k 33 

<x p / °C 


isotropic 

beam 

unit-cell 

analysis 

19.78x10® 

0 

5.35 

5.96x10® 

10x10® 

0 

beam 

theory 

19.78x10® 

0 

5.34 

5.77x10® 

10x10® 

0 

bimaterial 

beam 

unit-cell 

analysis 

74.29x10® 

30.20x1 0 3 

20.06 

8.47x10® 

30.73x10® 

-14.62x1 0‘ 3 

beam 

theory 

74.29x10® 

30.20x1 0 3 

20.06 

8.62x10® 

30.74x10® 

-14.63x1 O’ 3 

plain- 

weave 

textile 

beam 

unit-cell 

analysis 

27.76x10® 

0 

5.41 

9.21x10® 

12.66x10® 

-24.12x10® 

mosaic 

model 

71.48x10® 

0 

8.13 

8.14x10® 

4.39x10® 

0 


K 13> K 23 and ay are zero for all cases 


Table 9. Comparison of Beam CTE's for Plain-weave Textile Beam 



CTE’s from 
beam model 

CTE’s from 
continuum model 

% error 

a P x 10 6 °/C 

12.66 

11.30 

-10.46 

a M x 10‘ 6 °/C/m 

-24.12 

0 

100 

CLy °/C 

0 

0 

- 
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Table 10. Periodic Displacement BC's Imposed on the Lateral Faces of the Plate Unit-cell 


■ 


u(a,y)- 

u(0,y) 

v(a,y)- 

v (0,y) 

w(0,y) 

u(x,b)- 

u(x,0) 

v(x,b)- 

v(x,0) 

w(x,b)- 

w(x,0) 

AT 

1 . 

e M -1 
1 

a 

0 

0 

0 


0 

0 

2. 


0 

0 

0 

0 

b 

0 

0 

3. 

y xy0 M =i 

0 

a/2 

0 

b/2 

0 

0 

0 

D 

K X M =1 

az 

0 

-a 2 /2 

0 

0 

0 

0 

5. 

<=/ 

0 

0 

0 

0 

bz 

-b 2 /2 

0 

6. 

K xy M =l 

0 

az/2 

-ay/2 

bz/2 

0 

-bx/2 

0 

7. 

AT M =1 

0 

0 

0 

0 

0 

0 

1 
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Table 11. Non-zero [A], [fi] and [Z>] Coefficients for Example 1 (Isotropic Plate) using 
Finite Elements 



A„ 

xl0‘ 

A 12 

xlO 4 

A 22 

xlO 6 

Am 

xlO 6 

Du 

xlO’ 3 

d I2 

xlO 3 

D 2 2 

xlO' 3 

Du 

xlO' 3 

a/, a/ 
xl0‘/°C 

fiTEl- 20 
(FEA) 

2.810 

0.843 

2.810 

0.983 

15.320 

4.606 

15.320 

5.358 

10 

lamination 

theory 

2.810 

0.843 

2.810 

0.983 

15.310 

4.593 

15.310 

5.358 

10 


[A], [£] and [D] coefficients in SI units 


Table 12. Non-zero [A], [£] and [D] Coefficients for Example 2 (Bimaterial Plate) using 
Finite Elements 



An, A 22 
xlO 6 

A n 

xlO 6 

a 66 

xlO 6 

Bll> &22 
xlO 3 

b 12 

xlO 3 

B 66 

xlO 3 

uTE y-20 
(FEA) 

9.832 

2.043 

3.895 

- 0.563 

- 0.108 

- 0.228 

lamination theory 
for two plies 

9.832 

2.043 

3.895 

- 0.563 

- 0.108 

- 0.228 

lamination theory 
using 3-D elastic 
constants 

9.844 

2.048 

3.899 

0 

0 

0 



B a, D 22 

xlO 3 

D j2 

xlO 3 

d 66 

xlO 3 

a/, a/ 
xio 6 rc 

IV. P/ 
rc/m 

ftTEx- 20 
(FEA) 

53.590 

11.149 

21.220 

17.800 

0.170 

lamination theory for 
two plies 

53.573 

11.131 

21.220 

17.814 

0.170 

lamination theory 
using 3-D elastic 
constants 

53.762 

11.183 

21.293 

8.190 

0 


651 










































































Table 13. Non-zero [A], [B] and [D] Coefficients for Examples 3 and 4 (Unidirectional 
Composite) using Finite Elements 



Example 3 


fiTEx - 20 
(FEA) 

0.690 

0.149 

0.496 

0.177 

Halpin-Tsai eqns. 
and lamination 
theory 

0.673 

0.109 

0.363 

0.113 

fiTE %- 20 
(FEA) 

0.452 

0.062 

0.285 

0.114 

Halpin-Tsai eqns. 
and lamination 
theory 

0.444 

0.039 

0.151 

0.045 


Example 4 


a/x 

10 6 

a/x 

10 6 

15.489 

26.184 

15.625 

55.112 

7.378 

13.188 

6.774 

34.239 



Example 3 


Example 4 


Dji D 12 D 2 2 I D 66 


xlO 6 xlO 6 xlO 6 xlO 6 10 


fiTEx-20 

(FEA) 


Halpin-Tsai 
eqns. and 
lamination 
theory 



fiTEx- 20 
(FEA) 


Halpin-Tsai 
eqns. and 
lamination 
theory 



[A], [B] and [D] coefficients in SI units 
















































Table 14. Non-zero [A], [5] and [£>] Coefficients for Examples 5, 6 and 7 using Finite 
Elements 6 



Example 5 


Example 6 


ftTEx- 20 
(FEA) 


lamination theory 
using 3-D constants 


fiTE%- 20 
(FEA) 


lamination theory 
using 3-D constants 


Example 7 


ftTEx- 20 
(FEA) 


lamination theory 
using 3-D constants 



12.750 


15.040 


15.362 


16.020 


All 

xl0‘ 

xl0‘ 

B„ xlO 3 

0.598 

0.511 

0 

0.508 

0.483 

0 

3.332 

1.237 

0 

1.461 

1.237 

0 

1.065 

1.228 

0.515* 

0.756 

1.228 

0 



Example 6 


Example 7 


fiTEx- 20 
(FEA) 


lamination theory 
using 3-D constants 


In example 7, B }2 = -B„ and 


[A], [B] and [£>] coefficients in SI units 



653 





















































Table 15. SAM Results for Examples 1 and 2 



Example 2 pTEy^ 10 

(bimaterial (SAM) 

medium) ' 

exact 

solution 


E x , E y 
(GPa) 


36.03 



E t 

(GPa) 

G xz , G yz 
(GPa) 

G X y 

(GPa) 

v«, 

v*, 

8.72 

2.48 

15.23 

0.599 

0.183 

9.79 

2.48 

15.23 

0.312 

0.208 


Table 16. SAM Results for Examples 3 and 4 



ft rex-io 
(SAM) 

(uiU directional 


composite) 


Example 4 

(unidirectional 

composite) 


(FEA) 


rule of 
mixt./Halpin- 
Tsai eqns. 


ftTEx-10 

(SAM) 


fiTEx-20 

(FEA) 


rule of 
mixt./Halpin- 
Tsai eqns. 


e l 

(GPa) 

E T 

(GPa) 

Glt 

(GPa) 

64 

40.51 

11.17 

63.55 

36.48 

12.93 

64 

34.55 

11.26 

43.23 

22.11 

4.43 

43.12 

18.15 

5.59 

43.40 

14.79 

4.45 


Gj-j v LT Vjt 
(GPa) 


8.36 0.245 0.300 


9.94 0.300 0.233 


0.300 0.300 



3.06 0.159 0.237 




3.92 0.242 0.222 


0.260 0.252 










































Table 17. SAM Results for Examples 5, 6 and 7 




E x , E y 
(GPa) 

E z 

(GPa) 

G XZ> G yz 

(GPa) 

G xy 

(GPa) 

v«, 

v 

v xy 

Example 5 

(plain* 

weave) 

fiTEx - 10 
(SAM) 

11.52 

6.48 

1.60 

1.60 

0.396 

0.171 

fiTEx- 20 
(FEA) 

11.55 

6.26 

1.94 

2.12 

0.399 

0.186 

Dasgupta 

results 

14.38 

6.25 

1.94 

3.94 

0.463 

0.167 

Example 6 
(plain- 

weave) 

fiTEx - 10 
(SAM) 

64.49 

11.37 

3.97 

4.26 

0.398 

0.027 

f*TEx-2 0 
(FEA) 

58.27 

! 

10.92 

4.32 

4.84 

0.378 

0.098 

TEXCAD 

64.38 

11.49 

5.64 

4.87 

0.396 

0.027 

Example 7 
(5-harness 
weave) 

fiTE) c-10 
(SAM) 

65.61 

11.14 

3.66 

4.44 

0.363 

0.032 

fiTEx- 20 

(FEA) 

62.51 

11.02 

4.45 

4.80 

0.349 

0.047 

TEXCAD 

66.33 

11.51 

4.93 

4.89 

0.342 

0.034 
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(a) 


(b) 


Figure 2. (a) unidirectional composite; (b) unit-cell for the composite. 
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Figure 3. Yam pattern in a plain weave preform (unit-cell boundary in dotted lines). 
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(a) 


(b) 


Figure 5. (a) Layered medium; (b) bimaterial beam. 
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Figure 6. Textile beam: (a) undeformed unit-cell and deformation under: (b) unit extensional 
strain; (c) unit curvature; (d) unit shear strain, top and bottom surfaces are traction free; (e) unit 
shear strain, tractions allowed on top and bottom surfaces, (not to scale) 
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100 


80 - 

60 - 

40 - 



ta 

c 

= -40 - 


~ -60 - 
-80 - 
- 100 - 



v=o 

V^SxlO 5 

V^lxlO 6 



- 1.5 -1 - 0.5 0 0.5 1 


axial force (P) 


xlQ 3 


Figure 8. Beam failure envelope based on yarn failure. 



Figure 9. Beam failure envelope based on matrix failure. 
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bending moment (M) 


.5 

axial force (P) 


Figure 10 . Comparison of failure envelopes based on yarn failure. 



axial force (P) 


xlO 5 


Figure 11 . Comparison of failure envelopes based on matrix failure 





Figure 12. Hierarchy of discretization for a unit-cell to implement SAM: (a) unit-cell 
(macrolevel); (b) slice (mesolevel); (c) element (microlevel). 
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